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Introduction 

The  low  speed  performance,  safety  and  handling  qualities  of  a 
high  span  loading  aircraft  depend  on  the  structure  and  stability  of  the 
vortex  flow  created  by  the  aircraft.  In  a previous  work  by  Matol  (1975)^ , 
a preliminary  study  was  made  of  the  various  flow  field  components  and  of 
existing  models  of  the  wing-vortex  interactions  in  the  entire  flow.  This 
study  has  shown  that  existing  models  are  severely  limited  to  the  extent  to 
which  they  incorporate  the  effects  of  leading-edge  separation,  trai ling- 
edge  separation,  vortex  bursting,  and  complicated  flight  configurations. 
Therefore,  the  present  program  was  initiated  to  formulate  a three-dimensional 
lifting  surface  theory  for  the  steady,  symmetric  flow  about  a flat  plate 
delta  wing  at  moderate  angles  of  attack.  This  simple  problem  was  investigated 
to  ensure  the  tractabllity  of  the  problem  while  retaining  the  features  of 
flows  with  leading-edge  vortices.  The  study  was  restricted  to  the  near 
field  vortex-wing  interactions  although  the  techniques  employed  are  quite 
general . 

Historically,  lifting  surface  theories  have  been  developed  with 
several  underlying  assumptions.  First  one  assumes  that  the  flow  is  inviscid 
outside  of  the  immediate  neighborhood  of  the  wing-body  combination.  Secondly, 
one  generally  assumes  that  the  flow  is  incompressible.  The  second  restric- 
tion can  be  eliminated  by  scaling  for  the  linear,  steady  problem.  This 
solution  procedure  can  also  be  extended  to  solve  the  unsteady  problem. 

Finally,  traditional  lifting  surface  theories  are  linearized,  l.e.,  perturba 
tion  velocities  created  by  the  presence  of  the  wing  are  assumed  to  be  small 
compared  to  the  free  stream. 

In  a linearized  lifting  surface  theory,  the  lift  is  linear  as 


a function  of  the  angle  of  attack  for  the  small  angles  considered.  Also, 
the  linear  formulation  precludes  the  formation  of  free  vortex  sheets  except 
at  the  trailing  edge. 

These  are  reasonable  restrictions  for  high  aspect  ratio  wings 
at  low  angles  of  attack  employed  on  commercial  aircraft  and  military 
transports.  However,  for  fighters  with  highly  swept  wings  operating  at 
moderate  to  high  angles  of  attack,  this  model  of  a lifting  surface  is  no 
longer  applicable  for  the  entire  flight  regime. 

It  has  been  historically  observed  that  leading-edge  vortices 
appear  at  a moderate  angle  of  attack  of  approximately  10  degrees  for  sweep 
angles  of  approximately  60  degrees.  A description  by  Ornberg  (1954r  of 
such  a flow  is  given  in  Figure  1.  For  less  highly  swept  wings,  the  formula- 
tion of  the  leading-edge  vortices  is  delayed  to  higher  angles  of  attack. 

These  leading-edge  vortices  are  important  because  they  induce  high  velocities 
on  the  wing  and  control  surfaces  in  their  vicinity,  and  they  also  affect  the 
overall  flow  field  for  other  aircraft  in  their  vicinity. 

The  problem  with  leading-edge  vortices  has  aroused  considerable 
interest  with  the  advent  of  supersonic  aircraft  where  sweep  was  employed  to 
decrease  the  drag  at  supersonic  speeds.  An  early  attempt  to  study  this 
problem  was  made  by  Brown  and  Michael  (1955)^.  They  solved  the  problem  for 
a simplified  model  of  a flat  plate  delta  wing  at  moderate  angles  of  attack 
under  the  additional  restriction  of  slender  body  theory,  i.e.,  the  deriva- 
tives in  the  flow  direction  were  Ignored  and  only  those  In  the  transverse 
plane  were  considered.  Slender-body  theory  precludes  the  inclusion  of 
subsonic  traillng-edge  effects  and  Is  limited  to  the  study  of  simple 
geometries.  Brown  and  Michael  used  the  following  approximation  for  the 
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Figure  1.  Schematic  sketches  showing  flow  on  suction  side  of 
70°  flat  plate  delta  wii.3  at  a*15°  [after  tirnberg  (1954)  ] . 


leading- edge  vortex  sheet.  They  modeled  the  vortex  core  by  a line  vortex 
whose  strength  increased  linearly  along  its  axis,  which  was  consistent  with 
their  assumption  of  conical  flow.  The  vortex  was  fed  by  a cut,  i.e.,  a 
feeding  sheet  from  the  leading  edge  which  was  restricted  to  the  cross-flow 
plane.  This  model  of  the  vortex  sheet  will  be  referred  to  as  a vortex-cut 
combination. 

4 

Later,  Smith  (1966)  refined  the  Brown  and  Michael  model  of  the 
leading-edge  vortex  sheet  to  include  a representation  of  an  actual  force- 
free  vortex  sheet  as  well  as  the  vortex  core. 

However,  all  of  these  early  presentations  had  the  restrictive 

assumptions  of  conical,  slender  body  theory.  The  restriction  of  conical ity 
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has  been  removed  by  Smith  (1957)  and  others  for  several  geometries. 

Thus,  the  final  step  in  modeling  the  fully  three-dimensional  wing- 
body  configuration  was  to  remove  the  assumption  of  slender  body  theory. 

That  is,  a fully  three-dimensional  lifting  surface  theory  with  leading-edge 
vortices  would  have  to  be  formulated. 

To  extend  standard  linearized  lifting  surface  theories  to  include 
the  effects  of  leading-edge  vortices,  one  must  1)  allow  the  presence  of 
leading-edge  vortices  and  2)  include  non-linear  terms  in  the  pressure 
calculations.  These  two  changes  add  to  the  complexity  of  the  original 
problem  in  that  1)  the  location  of  the  leading-edge  vortex  sheets  is  origin 
ally  unknown  and  2)  for  a given  location  of  the  vortex  sheets,  the  require- 
ment of  zero  load  --  the  Kutta  condition  --  at  the  trailing  edge  becomes 
nonlinear  in  terms  of  the  velocity  or  vorticity  components. 

Linearized  lifting  surface  theories  have  been  developed  along 
two  distinct  lines.  First,  there  are  the  finite  element  methods  where  the 
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wing-body  combination  is  replaced  by  a finite  number  of  discrete  vortex 
elements  and  their  strengths  are  determined  by  satisfying  appropriate 
boundary  conditions.  The  alternate  method  is  to  employ  a finite  set  of 
loading  functions  whose  coefficients  are  then  chosen  to  satisfy  the  boundary 
conditions. 

Although  the  linearized  method  based  on  discrete  vortex  elements 
has  received  considerable  attention,  it  contains  several  inherent  disadvan- 
tages when  compared  with  the  alternate  approach  using  continuous  loading 
functions.  First,  there  is  an  infinite  discontinuity  in  the  downwash 
between  the  discrete  panels  or  at  the  vortex  element.  This  makes  the  solu- 
tion highly  sensitive  to  the  location  of  the  collocation  points  where  the 
velocity  is  being  matched.  Secondly,  finite  element  methods  have  tradition- 
ally had  stability  problems.  Artifices  have  been  introduced  to  ensure 
stability  in  given  instances,  but  such  modifications  introduce  further 
arbitrariness  into  the  final  solution.  Finally,  there  is  the  difficulty 
in  modeling  the  local  behavior  at  the  wing  edges  by  a choice  of  discrete 
elements. 

For  the  nonlinear  problem  with  leading-edge  vortices,  one  can 
also  attempt  to  model  the  wing  and  wake  by  finite  elements  or  by  continuous 
functions.  Due  to  the  difficulties  often  encountered  with  finite  element 
methods,  a three-dimensional  lifting  surface  theory  employing  continuous 
loading  functions  was  sought.  Previous  work  has  been  done  in  this  area  by 
Nangia  and  Hancock  (1968).**  Therefore,  the  model  of  Nangia  and  Hancock  was 
employed  as  a basis  for  this  investigation. 

Basically,  the  Nangia  and  Hancock  model  attempted  to  add  a Brown 
and  Michael  type  of  leading-edge  vortex  model  to  a conventional  lifting 
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surface  model.  Whereas  the  linearized  problem  can  be  solved  In  one  step 
as  the  solution  of  a system  of  simultaneous  linear  equations,  the  nonlinear 
problem  with  leading-edge  vortices  must  be  solved  by  an  iterative  procedure. 
The  method  that  Nangia  and  Hancock  employed  was  simple  in  theory  but  required 
extensive  computation.  It  was  hoped  that  with  the  advance  of  computational 
speed,  that  it  would  be  possible  to  reduce  the  original  lengthy  computations 
to  a reasonable  length.  Although  several  deficiencies  of  their  model  will 
be  noted  later,  their  model  was  used  as  a basis  to  take  advantage  of  the 
expertise  that  they  had  gained  in  their  computations.  Some  of  the  efforts  to 
refine  their  computations  will  be  discussed  later. 

The  basic  governing  equations  are  similar  to  those  employed  by 
Nangia  and  Hancock,  but  a new  iterative  procedure  has  been  developed.  This 
was  done  to  alleviate  several  shortcomings  in  the  original  procedure.  First, 
the  Nangia  and  Hancock  procedure  is  not  self-starting.  Secondly,  the  method 
requires  considerable  judgment  in  selecting  the  values  for  the  next  iteration. 

The  present  formulation  will  be  given  in  the  next  section. 

Problem  Formulation 

The  choice  of  axes  is  presented  in  Figure  2.  The  planform  is 
presently  considered  to  be  in  the  x-y  plane.  The  problem  can  also  be  treated 
if  the  lifting  surface  is  non-planar  if  a spanwlse  coordinate  is  used  to 
replace  y.  The  planform  can  be  completely  general.  If  there  is  no  side- 
slip, the  clean  configuration  can  often  be  considered  symmetric  and  the 
flow  field  can  be  completely  described  by  satisfying  the  boundary  conditions 
on  one  side  only.  The  boundary  conditions  on  the  other  side  are  automatically 
satisfied  by  symmetry.  To  facilitate  comparison  with  the  results  of  Nangia 
and  Hancock,  the  delta  wing  of  unit  aspect  ratio  was  chosen.  See  Figure  3 
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for  the  representation  of  the  wing,  leading-edge  vortices  and  wake. 

The  governing  equation  in  three-dimensional  invlscid,  irrotatlonal , 
steady  flow  about  a wing-body  combination  is  Laplace's  equation.  The  solu- 
tion can  be  formulated  as  an  integral  equation  over  the  boundary  of  the 
aircraft  configuration  and  regions  of  shed  vorticity.  There  are  several 
equivalent  formulations  for  the  solution,  but  vortex  sheets  are  used  to 
represent  the  wing  and  wake  in  this  formulation.  The  velocity  distribution 
can  then  be  given  in  the  following  vector  form 


where 


v(r)  = - 


1_ 

4?r 


S' 


Y X (f1  - r)  dS' 
!?•  - r|3 


0) 


r'  - r = (x ' - x)i  + (y'  - y)j  + (z*  - z)k 
Y = Yxi  + Yy  J + Yzk 

^ ^ A A 

v = ui  + vj  + wk 

S'  is  the  surface  of  integration,  y is  the  vorticity  vector,  v Is  the 
perturbation  vector  which  does  not  include  the  unit  free  stream,  and  r is 
the  radius  vector  from  the  origin.  The  velocities  are  nondl mens Iona 11  zed 
with  respect  to  the  free  stream  and  the  distances  are  nondlmenslonallzed 
with  respect  to  the  root  chord. 

Since  the  vorticity  lies  In  the  plane  of  the  wing  and  wake,  the 
vorticity  representing  the  wing  consists  of  only  two  nonzero  components, 

Y and  y . Continuity  of  vorticity  can  be  represented  by 

a y 


3Yy  x ^x 

"57"  3x 


(2) 
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t. 


This  implies  that  only  one  component  of  vorticity  need  be  given  on  the 
wing  and  that  the  other  component  can  then  be  calculated  from  Equation  2. 

In  the  present  formulation,  the  vorticity  on  the  wing  is  divided  into  two 
parts.  First,  there  is  a portion — subscript  1 — which  behaves  like  the 
traditional  bound  vortices  and  does  not  feed  into  the  leading-edge  vortex 
sheets.  Secondly,  there  is  a portion— subscript  2— which  feeds  the  leading- 
edge  vortices. 

Yy  s Y = Yi  + (3) 

YX  5 6 * 61  + 52 

Y and  6 are  used  to  represent  the  vorticity  components  to  simplify  comparison 
with  the  original  Nangia  and  Hancock  report.  The  contributions  are  chosen 
so  that  y.|  and  <5-|  fall  to  zero  at  the  leading  edge,  while  Y£  and  are 

chosen  so  that  the  vorticity  is  perpendicular  to  the  leading  edge.  This  is 
necessitated  by  the  Brown  and  Michael  model  employing  a vortex-cut  combination 
to  ensure  finite  velocities  at  the  leading  edge.  A more  general  model  employ- 
ing the  leading-edge  vortex  sheet  formulated  by  Smith  would  utilize  a 
vorticity  distribution  which  resulted  in  a no-load  condition  at  the  leading 


The  functional  forms  of  the  wing  vorticity  are  modeled  after  those 


employed  by  Nangia  and  Hancock. 


n m 


6,  = l £ a xq  (y/kx) ( 1 - (y/kxT) 
1 p*0  q*0  p,q 

x S fx2+y2Kq-1)/2 


2p+l 
2x  c 


*2  ■ OT7?  ji  q K 


(4) 


while  the  leading-edge  vortex  strength  Is  defined  by 


T s [ 9 xq  (5) 

q=l  H 

See  Figure  4 for  a representation  of  bound  vortlclty  62*  The  coefficients 

a and  gn  are  the  unknown  coefficients  of  the  vortlclty  functions,  k Is 
p.q  q 

the  semispan  at  the  trailing  edge.  The  functions  have  been  chosen  to  satisfy 
the  aforementioned  restrictions  and  the  obvious  symmetry  constraints. 

The  wake  has  been  considered  to  be  flat  and  an  extension  of  the 
vorticlty  strength  at  the  trailing  edge  as  In  linear  lifting  surface  theory. 
Consequently,  the  trailing  wake  adds  no  new  unknown  parameters. 

Finally,  the  location  of  the  leading-edge  vortex  is  defined  by 
the  polynomials 


y»  ■ J,  v x 


*v  ■ J,  ■» x 


(6) 


This  Introduces  the  final  unknowns  gyy,  gzy 

If  there  were  no  numerical  difficulties  associated  with  this  prob- 
lem, one  could  obtain  Improved  accuracy  by  Increasing  the  number  of  modes 
and,  consequently,  increasing  the  number  of  unknowns.  However,  for  the 
given  set  where  both  the  leading-edge  vortex  strength  and  vortex  locations 
are  defined  by  polynomials,  a large  number  of  modes  also  Increases  the 
problems  related  with  the  numerical  analysis  necessary  to  solve  the  result- 
ing equations. 

After  the  mode  shapes  have  been  defined.  It  is  necessary  to  deter- 
mine appropriate  governing  equations  to  determine  the  unknown  coefficients. 
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Figure  4.  Representation  of  bound  vorticity 
feeding  leading-edge  vortex. 


The  governing  equations  are  the  no-flow  condition  through  the  wing  and  the 
no-load  condition  on  the  free  vortex  sheet  and  on  the  leading-edge  vortex- 
cut  combination. 

The  downwash  condition  becomes 


w * - sina  (7) 

where  a is  the  angle  of  attack  and  w is  the  downwash  Induced  by  the  vortlclty 
distribution.  This  requires  the  evaluation  of  the  w component  of  the  Integral 
In  Equation  1 at  a set  of  collocation  points.  The  sinusoidal  distribution 
of  collocation  points  of  Multhopp  (1950)^  modified  for  the  delta  wing  Is 
given  In  Figure  5.  A large  amount  of  computational  time  is  consumed  by  the 
calculation  of  the  contribution  from  the  wing  surface  since  the  denominator 
contains  a singularity  at  the  collocation  points.  Figure  6 gives  the 
Integration  regions  employed  for  the  surface  Integration.  Each  region  was 
covered  by  a 24x24-po1nt  Gaussian  quadrature.  Fortunately,  this  computation 
does  not  require  any  Iterations. 

The  no-load  condition  on  the  trailing  sheet  Is  satisfied  approxi- 
mately. First,  the  Kutta  condition  at  the  trailing  edge  requires  the 
vortlclty  vector  to  be  parallel  to  the  velocity  vector. 


Y(l.y) 

= tanB 

6(1. y) 


v 

cosa  + u 


(8) 


This  equation  Is  employed  at  collocation  points  along  the  trailing  edge. 


Secondly,  aft  of  the  trailing  edge,  the  wake  Is  considered  to  be  the  same 


as  In  the  linearized  lifting  surface  theory,  since  the  major  cause  of 
nonlinearity  Is  the  presence  of  the  leading-edge  vortices  which  Induce 
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Figure  5.  Wing  collocation  points. 
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high  spanwise  velocities  on  the  planform. 

Finally,  the  condition  of  no-load  on  the  vortex  cut  is  formulated 
on  the  right-hand  vortex  as  an  extension  of  the  Brown  and  Michael  model. 

The  force  components  in  the  y and  z directions  are  given  by  and  Fz> 
respectively. 

dy  (9) 

where  w1 , v1  are  the  velocities  induced  by  the  vorticity  distributions 
excluding  the  contributions  of  the  right-hand  vortex  on  itself.  The  calcula- 
tions of  the  velocities  v^  and  w1  at  collocation  points  along  the  leading- 
edge  vortex  again  require  the  evaluation  of  the  surface  integral  in  Equation 
1.  The  integrand  is  no  longer  singular  and  two  24x24-point  Gaussian  quadra- 
tures are  applied  over  the  wing  surface. 

The  procedure  to  calculate  the  unknowns  Is  now  described.  First, 
the  number  of  modes  to  be  employed  must  be  specified.  Once  the  number  of 
modes  and  consequently,  the  number  of  unknown  vorticity  coefficients  have 
been  determined,  it  is  necessary  to  choose  an  Iterative  scheme  to  satisfy 
the  system  of  equations  provided  by  the  downwash  condition,  the  Kutta 
condition  and  the  no-force  condition.  For  a given  vortex  location,  the 
downwash  condition  is  linear  In  terms  of  the  vorticity  distribution  while 
the  equations  provided  by  the  Kutta  condition  and  the  no-force  condition  are 
nonlinear  In  all  parameters. 

A preliminary  formulation  followed  the  procedure  of  Nangla  and 
Hancock  and  only  attempted  to  reduce  the  amount  of  computation  Involved 
during  each  step.  Their  method  can  be  summarized  as  follows: 
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1)  Evaluate  the  contributions  to  the  downwash  from  the  wing 
vorticity. 

2)  Specify  the  initial  position  of  the  vortex.  For  example,  this 
can  be  furnished  by  the  Brown  and  Michael  model. 

3)  Specify  an  initial  distribution  of  B(y)  such  that  it  does 
not  exceed  tan'^(y). 

4)  Evaluate  the  downwash  integral  due  to  the  vortices  and  wake. 

5)  Formulate  the  Kutta  condition  as  linear  equations  in  the 
unknown  vorticity  distribution  at  the  collocation  points 
using  the  specified  tanB. 

6)  Solve  the  simultaneous  equations  from  the  downwash  and  Kutta 
conditions  for  the  vorticity  coefficients  ap  and  gq. 

7)  Use  the  results  of  the  previous  steps  to  determine  new  values 
for  B. 

8)  If  the  new  value  of  8 is  greater  than  or  equal  to  the  Initial 
value  in  step  3),  go  on  to  the  next  step;  for  smaller  values 
of  B,  iterate  between  4)  and  7)  until  the  value  of  8 converges. 
This  conditional  form  for  calculating  B is  employed,  because 
Nangla  and  Hancock  found  that  the  procedure  otherwise  furnished 
unrealistic  results. 

9)  Calculate  the  forces  on  the  right-hand  vortex  using  the  assumed 
vortex  position  and  the  previously  obtained  vorticity  distribu- 
tions. If  the  force  exceeds  a prescribed  tolerance,  then  the 
slopes  of  the  new  vortex  location  are  calculated  from 


19 


110) 


where  d is  chosen  small  enough  to  prevent  divergence  of  the 
method,  e.g.,  d = .01.  Since  the  forces  have  been  normalized, 
the  vortex  movement  is  restricted  to  d. 

10)  Iterate  between  steps  4)  and  9)  until  the  forces  become  suffi- 
ciently small.  The  solution  has  ti.an  been  obtained  for  the 
given  flow  conditions. 

This  original  procedure  has  been  modified  in  the  present  procedure 
so  that  the  program  can  be  run  automatically  with  a minimum  of  stops  to 
input  new  data. 

Instead  of  specifying  6 in  Step  3,  the  downwash  condition  was 
linearized  with  respect  to  the  vorticity  coefficients  and  the  initial  approx- 
imation for  the  coefficients  an  „ and  g are  found  directly  from  the  downwash 

p.q 

equations.  Then  $ can  be  calculated  by  employing  Equation  8 

The  calculation  of  the  forces  Fy  and  Fz  involve  integrations  over 
the  wing  surface  to  obtain  the  contribution  to  the  v and  w velocity  components 
at  the  collocation  points  along  the  leading-edge  vortex.  Such  integrations 
are  time-consuming  in  the  present  scheme.  Since  a large  number  of  iterations 
are  required,  a substantial  savings  could  be  gained  by  developing  an  iterative 
scheme  with  faster  convergence. 

In  the  two-dimensional  theory  of  Brown  and  Michael,  one  encounters 
a similar  problem  in  solving  for  the  vortex  location,  yy  and  zy,  in  the  cross- 
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flow  plane  by  satisfying  the  no-force  condition.  Brown  and  Michael  originally 
solved  the  problem  indirectly  by  assuming  values  of  the  vortex  location  and 
satisfied  the  no-force  condition  by  trial  and  error.  Later,  Pullin  (1973) 
developed  a Newton's  method  for  the  Smith  type  model,  which  included  the 
Brown  and  Michael  solution  as  a degenerate  case.  Trial  runs  of  Pul  Tin's 
program  indicated  that  the  Newton's  method  "converged"  in  approximately  four 
iterations  for  the  Brown  and  Michael  model. 

A Newton's  method  has  several  advantages  over  the  Nangia  and  Hancock 
procedure.  First,  the  scheme  is  amenable  to  automatic  iteration  without 
operator  interference.  Secondly,  Newton's  method  has  been  shown  to  be 
successful  for  a wide  variety  of  problems.  The  Nangia  and  Hancock  procedure, 
on  the  other  hand,  has  neither  of  these  virtues.  First,  from  Equation  10, 
the  extent  of  the  vortex  movement  is  always  ' d ' , which  must  be  specified  by 
the  operator.  Secondly,  by  comparing  Equations  9 and  10,  one  sees  that  the 
Nangia  and  Hancock  procedure  for  moving  the  vortex  is  equivalent  to  a Newton's 
method  in  which  the  velocities,  Vj  and  Wj , are  considered  constant  with 
respect  to  the  vortex  movement.  The  results  of  numerical  experiments  dis- 
cussed in  the  next  section  indicate  that  the  ignored  derivatives  are  impor- 
tant. Consequently,  although  the  Nangia  and  Hancock  method  appears  to  work 
for  the  simple  case  they  considered,  the  procedure  is  not  generally  valid 
and  has  a convergence  rate  which  depends  on  the  applicability  of  the  proce- 
dure and  the  judgment  of  the  operator  In  selecting  values  for  ' d ' . Therefore, 
a Newton's  method  was  developed  to  locate  the  new  vortex  position  from  the 
forces  calculated  In  the  preceding  Iteration.  This  replaces  step  9)  in  the 
Nangia  and  Hancock  procedure. 
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Results  of  the  Program 

The  problem  was  formulated  for  the  unit  aspect  ratio  delta  wing 
at  an  angle  of  attack  of  14.3°.  Four  spanwlse  modes  and  5 chordwise 
modes  were  employed  to  describe  the  surface  vorticity,  since  these  were  the 
number  of  nodes  employed  by  Nangia  and  Hancock.  This  choice  for  the  vorticity 
representation  corresponds  to  twenty  unknown  a *s  and  five  g 's  In  Equation 
4 for  a total  of  twenty-five  unknown  vorticity  coefficients.  The  Initial 
location  of  the  vortex  was  obtained  from  the  conical,  slender  model  of 
Brown  and  Michael.  An  approximation  for  the  vorticity  coefficients  was  then 
obtained  by  satisfying  the  downwash  condition  at  twenty-five  collocation 
points.  These  collocation  points  included  the  twenty  interior  points  In 
Figure  5,  as  well  as  an  intermediate  row  of  five  points  at  x * .86.  This 
initial  solution  serves  three  primary  purposes.  First,  it  can  be  compared 
with  the  results  of  Brown  and  Michael  and  similar  slender  body  theories 
to  check  the  convergence  of  the  method  as  a function  of  the  number  of  modes. 
Secondly,  it  can  provide  insight  on  better  choices  for  the  vorticity  functions. 
For  example,  in  linearized  lifting  surface  theory,  the  modes  are  chosen  to 
approximate  the  two-dimensional  flat  plate  results  and  the  elliptic  wing 
result.  Unfortunately,  It  Is  more  difficult  to  obtain  such  a simple  choice 
of  vorticity  functions  for  a leading-edge  vortex  problem  as  the  suction  peak 
location  is  a function  of  the  vortex  core  location,  which  Is  originally  un- 
known. Hopefully,  a better  choice  of  vorticity  functions  can  be  developed 
in  the  future  from  the  results  of  this  Investigation.  Finally,  this  Initial 
set  of  vorticity  coefficients  is  required  to  calculate  an  Initial  choice  of 
3,  which  Is  then  used  to  formulate  the  quasl-1 inear  form  of  the  Kutta 
condition  at  the  trailing  edge  from  Equation  8. 


The  original  results  for  the  wing  loading  versus  non-dimensional i zed 
span  are  plotted  in  Figure  7 for  two  chordwije  stations.  These  results 
are  also  compared  with  those  of  th*  Brown  and  Michael  model  for  the  equiva- 
lent aspect  ratio  and  angle  of  attack.  The  present  calculated  values  are 
similar  to  the  slender  body  results  but  the  loading  decreases  towards  the 
trailing  edge  as  a result  of  the  finite  planform. 

Secondly,  the  Kutta  condition  was  applied  at  the  five  collocation 
points  at  the  trailing  edge  in  Figure  5.  The  downwash  condition  previously 
applied  at  the  twenty  interior  collocation  points  was  used  with  the  five  new 
equations  to  provide  enough  constraints  to  calculate  the  new  set  of  twenty- 
five  unknown  coefficients.  The  resulting  loading  distributions  are  shown  in 
Figure  8.  Note  that  the  loading  goes  to  zero  at  the  trailing  edge  in  this 
moiel  which  is  not  possible  in  slender  body  models.  Also,  note  the  peaks 
in  the  loading  which  are  not  possible  in  linearized  models.  Some  features 
of  the  model  can  be  observed  by  examining  a plot  of  the  leading-edge  vortex 

HP 

strength  for  the  same  case  in  Figure  9,  Strictly  speaking,  should  be 
zero  at  the  trailing  edge  in  this  model,  since  the  leading  edge  vortex  Is 
only  fed  along  the  leading  edge  of  the  planform.  This  condition  was  included 
In  a later  run,  but  since  it  only  affected  the  shape  of  the  curve  near  the 
trailing  edge,  this  constraint  was  not  incorporated  In  this  early  model. 
However,  in  later  models,  when  the  force  Is  to  be  calculated  at  points  near 
the  trailing  edge,  this  imposed  boundary  condition  becomes  Important  and 
would  be  best  handled  by  a choice  of  modes  which  Incorporated  this  condition 
automatically.  Although  the  magnitude  of  the  leading-edge  vortex  strength 
changed  only  slightly  with  the  application  of  the  Kutta  condition,  the  func- 
tional shape  became  more  pronounced.  Numerical  experiments  were  conducted 
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Figure  7.  Loading  calculated  by  satisfying  downwash  condition 
only.  Symbols  indicate  points  where  pressure  was  calculated. 

25  combined  modes. 
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Figure  9.  Leading-edge  vortex  strength.  Iterations  1 and- 2 
refer  to  number  of  applications  of  the  Kutta  condition. 

25  combined  modes. 


to  demonstrate  that  the  curvature  of  the  function  representing  the  vortex 
strenqth  was  not  a numerical  error  or  a result  of  the  collocation  points. 
Rather,  it  is  a feature  of  the  present  formulation  which  Is  alleviated  by 
further  iterations. 

Preliminary  results  indicated  large  values  with  alternating  signs 
for  the  unknown  vorticity  coefficients,  which  suggest  several  difficulties. 
First,  there  is  the  loss  of  significant  figures  In  the  final  answer  when 
two  large  numbers  are  differenced  to  obtain  a small  number.  Secondly,  they 
indicate  that  the  original  choice  of  vorticity  functions  was  probably  not 
a good  one.  Finally,  their  values  can  change  radically  from  one  iteration 
to  the  next.  This  is  important,  since  one  could  save  considerable  computa- 
tional effort  if  a quasi-linearization  could  be  employed  about  some  initial 
value.  For  example,  it  may  be  possible  to  Ignore  some  of  the  derivatives 
in  the  Jacobian  or,  for  slow  variations,  some  nonzero  terms  may  only  have 
to  be  calculated  once  for  several  iterations. 

The  present  set  of  equations  for  the  unknown  vorticity  coefficients 
is  also  sensitive  to  small  perturbations  In  the  boundary  conditions.  This 
is  due  to  the  fact  that  the  relatively  large  matrix,  25x25,  Is  not  diagonally 
dominant,  which  causes  numerical  difficulty  In  the  solution  procedure. 

To  reduce  these  numerical  difficulties  and  decrease  the  computation- 
al effort.  It  would  be  advantageous  to  utilize  the  minimum  number  of  modes 
which  would  adequately  describe  the  spanwlse  loading  distribution.  Therefore, 
the  computer  program  developed  to  calculate  the  vorticity  coefficients  from 
the  downwash  condition  was  used  to  approximately  determine  the  minimum  number 
of  spanwlse  modes. 

This  was  done  by  using  a quasi -conical  model  where  only  the  linear 
chordwise  modes  were  retained  in  the  set  of  vorticity  functions.  The  location 
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of  the  leading-edge  vortex  was  taken  from  the  Brown  and  Michael  model  and 
should,  therefore,  satisfy  the  no-force  condition.  This  model  is  quasi- 
conical,  in  the  sense  that  the  integration  is  terminated  at  the  trailing  edge 
of  this  subsonic  model  and  is  not  extended  to  infinity.  The  results  of  these 
calculations  should  only  match  the  Brown  and  Michael  solutions  approximately. 
The  pressure  distribution  obtained  using  three  and  four  spanwise  modes  and 
twenty  Interior  collocation  points  in  a least-square  method  is  compared  with 
the  conical  model  in  Figure  10.  The  leading-edge  vortex  strength,  which  is 
now  linear,  is  compared  in  Figure  11.  These  results  indicate  that  four  span- 
wise  modes  are  sufficient  to  adequately  describe  the  spanwise  loading  distribu- 
tion with  the  given  choice  of  loading  functions. 

A similar  study  was  not  conducted  to  determine  the  minimum  number 
of  chordwlse  modes,  because  the  original  modes,  given  In  Equation  4,  behave 
like  simple  polynomials  in  chordwise  direction.  Such  polynomials  are  asso- 
ciated with  loss  of  accuracy,  and  attention  was  focused  on  attempting  to 
develop  a better  chordwlse  dependence  rather  than  determining  their  minimum. 
However,  if  these  original  chordwlse  modes  were  to  be  used,  a least  squares 
method  should  be  employed  to  reduce  this  problem. 

Finally,  a program  was  developed  to  calculate  the  forces  on  the 
vortex  and  to  move  the  vortex  according  to  Equation  10.  This  program 
required  approximately  one  minute  to  calculate  the  forces  at  five  chordwlse 
stations.  To  prevent  divergence,  the  movement  had  to  be  small  which  implied 
a large  number  of  iterations.  Therefore,  a procedure  more  amenable  to  auto- 
matic iteration  was  sought. 

Thus,  a numerical  experiment  on  the  applicability  of  Newton's 
method  was  conducted  for  the  two-dimensional  model  of  Brown  and  Michael,  since 
it  was  felt  that  useful  information  could  be  obtained  on  the  force  Jacobian 
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Figure  11.  Leading-edge  vortex  strength 
for  3 and  4 spanwise  modes. 
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more  economically  in  two  dimensions  than  in  the  fully  three-dimensional  case. 
The  numerical  experiment  was  conducted  on  a slender  delta  wing  of  the  same 
aspect  ratio  and  angle  of  attack  as  in  the  present  case  for  different  leading- 
edge  vortex  strengths  and  locations  which  correspond  to  the  unknowns  in  the 
present  case.  Equation  9 was  employed  to  calculate  the  forces,  Fy  and  Fz. 

In  the  conical  Brown  and  Michael  model,  both  the  vortex  location,  yy  and  zy, 
and  the  vortex  strength,  r,  are  linear  in  the  chordwise  variable. 

Results  of  this  experiment  for  r(l.)  = .31  are  plotted  in  Figure 
12  for  the  spanwise  force,  Fy,  and  in  Figure  13  for  the  vertical  force,  Fz- 
The  forces  are  plotted  versus  vortex  location,  yy  and  zy,  at  x * 1.  The 
triangle  symbol  represents  the  point  where  the  lines,  Fy  = 0 and  Fz  = 0, 
intersect  to  define  the  stable  location  for  this  flow  condition.  The  two- 
dimensional  results  indicate  that  F?  is  monotonic  as  a function  of  both  yy 
and  zy.  Fy,  on  the  other  hand,  is  monotonic  over  much  of  the  plane,  but 
is  poorly  behaved  near  the  leading-edge.  This  did  not  preclude  the  use  of 
Newton's  method  on  the  Brown  and  Michael  model,  but  may  prove  to  be  a 
problem  in  the  three-dimensional  case,  since  Newton's  method  is  based  on  a 
linear  approximation.  However,  it  should  be  possible  to  minimize  this 
problem  by  avoiding  this  region  or  by  suppressing  some  of  the  derivatives 
or  by  limiting  changes  of  the  leading-edge  vortex  location  to  one  axis  at 
a time. 

Therefore,  a Newton's  procedure  was  developed  for  the  three- 
dimensional  case  to  calculate  the  new  location  of  the  vortex  based  on  the 
forces  calculated  In  the  preceding  Iteration.  This  modification  Improved 
the  rate  of  convergence  for  the  reduction  of  the  forces  for  a given  vorticity 
distribution.  However,  when  the  given  vorticity  distribution  was  altered 
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by  Iterating  on  the  downwash  and  Kutta  conditions,  the  large  changes  in  the 
vorticity  coefficients  neutralized  much  of  the  gain  in  the  convergence  rate. 

A large  number  of  Iterations  were  therefore  not  conducted,  because  a more 
economical  computing  procedure  must  be  developed  before  this  method  can  be 
used  economically.  As  is,  the  program  will  produce  answers  for  the  delta 
wing  problem  given  an  initial  vortex  location  and  outputs  from  supplementary 
programs  which  calculate  the  downwash  contribution  from  the  wing  and  an 
initial  approximation  for  the  leading-edge  vortex  strength. 

Conclusion 

In  conclusion,  further  effort  should  be  devoted  to  model  the 
leading-edge  vortex  problem  for  arbitrary  planforms  by  a continuous  vorticity 
model  with  unknown  coefficients,  since  the  role  of  vortex-wing  interactions 
becomes  Increasingly  crucial  for  high  technology  aircraft.  The  present 
procedure  can  be  easily  generalized  to  arbitrary  planforms  and  Include 
arbitrary  sources  of  free  vortices,  but  Is  restricted  by  the  computational 
effort. 

Due  to  the  advance  of  computational  power  and  some  minor  refinements 
employed  In  these  calculations,  the  calculation  of  the  downwash  contribution 
from  the  wing  has  been  reduced  from  the  two  and  a half  hours  on  an  Atlas® 
to  twelve  minutes  for  twenty  collection  points  and  25  unknown  vorticity 
coefficients  on  an  IBM  370-168.  This  downwash  calculation  need  only  be 
performed  once  for  a given  planform  and  a given  set  of  collocation  points, 
and  solutions  can  then  be  obtained  for  all  angles  of  attack  and  for  different 
models  of  the  trailing  wake,  leading-edge  vortex,  and  leading-edge  vortex 
sheet.  Thus,  It  was  originally  considered  justified  to  devote  this  much 
effort  to  this  downwash  calculation,  since  it  was  implied  by  Nangia  and 
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Hancock  that  this  calculation  would  require  the  major  portion  of  computational 
effort.  However,  the  present  investigation  indicates  that  the  program  for 
reducing  the  forces  while  satisfying  the  downwash  and  Kutta  Conditions  consumes 
an  equivalent  or  even  greater  effort.  Therefore,  an  important  conclusion  is 
that  several  important  modifications  must  be  utilized  to  make  the  present 
program  generally  useful.  Primarily,  it  would  be  desirable  to  reduce  the 
order  of  the  matrices  in  terms  of  the  numerical  instabilities  and  to  further 
reduce  the  computational  effort  to  improve  the  versatility  of  the  method. 

This  would  involve  improved  choices  for  the  vorticity  functions,  the  use  of 
a least  squares  method,  and  an  alternate  formulation  of  the  integral  equation. 

Presently,  four  spanwise  modes  of  the  Nangia  and  Hancock  type 
adequately  describe  the  spanwise  loading.  However,  due  to  the  large  values 
of  the  coefficients  with  alternating  signs,  better  accuracy  can  probably  be 
obtained  by  choosing  more  orthogonal  functions,  which  retain  the  square  root 
behavior  near  the  leading  edge.  On  the  other  hand,  the  chordwlse  variation 
must  be  altered  to  satisfy  the  dr/dx  = 0 condition  at  the  trailing  edge  and 
to  eliminate  the  numerical  difficulties  associated  with  simple  polynomials. 

The  least-squares  method,  mentioned  in  connection  with  determining 
the  minimum  number  of  spanwise  modes,  is  an  important  technique  for  alleviat- 
ing some  of  the  numerical  difficulties  associated  with  large  coefficients  of 
alternating  sign  and  with  large  matrices.  The  least-squares  technique  does 
not  attempt  to  solve  the  set  of  simultaneous  equations  at  a finite  number  of 
collocation  points  exactly.  Instead,  the  number  of  collocation  points  is 
chosen  to  be  greater  than  the  number  of  unknowns,  and  the  problem  is  formulated 
to  minimize  the  error  at  all  points.  This  type  of  formulation  results  in 
smoother  solutions,  i.e.,  smaller  values  of  coefficients,  and  tends  to  mini- 
mize local  anomalies.  This  method  also  allows  a smaller  number  of  modes  to 
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be  utilized  to  adequately  predict  the  flow  characteristics.  Finally,  the 
least-squares  procedure  results  In  a more  diagonally  dominant  matrix  than 
can  be  obtained  by  simply  applying  the  boundary  conditions. 

Finally,  an  alternate  formulation  of  the  Integral  equation  Is 
suggested.  Basically,  the  vortlcity  which  does  not  feed  the  leading-edge 
vortex  can  be  represented  by  a vortlcity  distribution  borrowed  from  a linear 
lifting  surface  theory.  Such  Integrations  In  lifting  surface  theory  routine- 
ly require  less  than  a minute,  and  this  alteration  would  dractlcally  reduce 
the  amount  of  computational  effort  required  to  calculate  the  velocity 
contributions  from  the  surface  vortlcity. 

These  modifications  have  been  suggested  by  the  results  of  the 
present  Investigation,  and  their  successful  Implementation  would  greatly 
enhance  the  versatility  of  a three-dimensional  lifting  surface  theory  with 
leading-edge  vortices. 
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SYMBOLS 


vorticlty  coefficient  related  to  surface  vorticity  only 
chord;  all  lengths  are  nond i mens ionali zed  by  chord 
scale  factor  to  limit  movement  of  vortex 
force  on  right-hand  vortex 

vorticity  coefficient  related  to  leading-edge  vortex 
unit  vector  in  x-direction 
unit  vector  in  y-directlon 

unit  vector  in  z-directlon;  semispan  at  trailing  edge 
dummy  index 
dummy  index 
dummy  index 

radius  vector  from  origin 
surface  of  integration 

free  stream  velocity;  all  velocities  are  nond 1 mens Iona 11 zed  by  U 
perturbation  velocity  component  In  x-dlrection 
perturbation  velocity  in  y-dlrectlon;  vector  v Is  total 
perturbation  velocity 

perturbation  velocity  component  In  z-dlrectlon 
chordwise  coordinate 

spanwlse  coordinate;  with  subscript  v,  represents  leading-edge 
vortex  location 

yertlcal  coordinate;  with  subscript  v,  represents  leading-edge 
vortex  location 
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I SYMBOLS  (continued) 

a angle  of  attack 

B shedding  angle  at  trailing  edge 

Y spanwlse  vortlclty  component;  vector  y Is  total  vortlclty 

r leading-edge  vortex  strength 

<5  chordwlse  vortlclty  component 
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